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LONG-TERM  GOALS 

To  understand,  characterize,  and  predict  lateral  and  vertical,  naturally-occurring  heterogeneities  that 
may  produce  discrete  acoustic  returns  at  low  grazing  angles  (i.e.,  "geologic  clutter")  in  a  mid-outer 
shelf  test  site  off  the  U.S.  (New  Jersey),  and  then  2)  to  conduct  precise  acoustic  reverberation 
experiments  at  this  site  to  understand,  characterize,  and  potentially  mitigate  the  geologic  clutter,  so  that 
the  false  alarms,  or  detects,  of  tactical  sonar  systems  encountered  in  this  marine  geologic  environment 
around  the  world  can  be  characterized  properly. 

OBJECTIVES 

The  primary  objective  of  this  effort  is  to  conduct  geostatistical  modeling  that  will  be  essential  for 
transitioning  from  STRATAFORM  products  to  Geoclutter  needs.  One  of  my  primary  STRATAFORM 
tasks  was  to  develop  a  means  to  realistically  interpolate  stratigraphic  architecture  from  limited  data, 
which  resulted  in  the  “SimStrat”  algorithm  (Goff,  2000).  Ultimately  this  work  is  intended  to  facilitate 
acoustic  reverberation  modeling  in  the  shallow  water  enviromnent.  With  the  Geoclutter  initiative,  a 
collaborative  effort  among  ONR  geophysicists  and  acousticians  to  understand  signal-like  reverberation 
in  the  natural  environment,  this  goal  is  on  the  verge  of  being  fully  realized.  The  New  Jersey 
STRATAFORM  natural  laboratory  has  also  been  chosen  for  the  Geoclutter  field  work,  in  part  because 
of  the  abundance  of  data  already  available,  and  also  because  of  the  complexity  of  the  seafloor  and 
subseafloor  environment.  In  particular,  numerous  buried  channels  just  meters  below  the  seafloor  are 
probable  candidates  for  “geoclutter”-type  returns.  This  “channels”  horizon  is  very  complex  in  nature, 
and  not  readily  modeled  within  the  SimStrat  framework.  The  goal  here  is  therefore  to  expand  my 
STRATAFORM  modeling  efforts  to  incorporate  channels-type  horizons. 

APPROACH 

SimStrat  is  intended  to  work  on  continuous  or  nearly  continuous  sets  of  stratigraphic  surfaces  (simple 
truncations  and  faults  could  be  handled  with  some  straightforward  adaptations).  While  this  is  true  of 
most  stratigraphic  horizons,  channeled  horizons  are  notable  and  critical  exceptions.  As  revealed  on  the 
New  Jersey  margin  seismic  data  (Figure  1),  the  channels  horizon  is  highly  discontinuous,  existing  only 
where  the  channels  themselves  exist.  Furthermore,  these  channels  clearly  form  a  dendritic  pattern  - 
which  is  not  reproducible  with  the  usual  suite  of  geostatistical  simulation  methods,  one  of  which  (the 
Fourier  method)  is  used  in  SimStrat.  A  different  approach  to  simulating  this  type  of  surface  is 
required. 
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A  channels-type  horizon,  with  dendritic  pattern  and  patchy  existence,  is  not  amendable  to  the  SimStrat 
algorithm  because  it  cannot  unconditionally  or  conditionally  be  simulated  by  conventional 
geostatistical  techniques.  An  entirely  new  method  for  interpolation  must  be  found  for  such 
morphologies.  Process-based  approaches  could  potentially  work,  but  the  goal  here,  as  with  SimStrat, 
is  to  base  the  conditional  simulation  only  on  geometric  and  statistical  considerations;  we  seek  speed 
and  simplicity  as  well  as  realism. 

In  the  New  Jersey  channels  horizon,  each  channel  can  be  considered  an  independent  surface  until  it 
merges  with  another  channel.  Owing  largely  to  this  observation,  I  believe  the  best  approach  to 
interpolating  a  channels  horizon  from  limited  sampling  (Figure  2  is  an  excellent  example)  is  to 
conditionally  simulate  each  channel  separately.  This  approach  will  greatly  simplify  the  problem  by 
breaking  it  down  into  manageable  parts,  but  it  will  also  require  sufficient  sampling  density  so  that  each 
primary  channel  can  be  individually  discerned  and  traced  through  the  interpolating  region.  On  the 
other  hand,  with  any  less  data  an  interpolation  of  channel  morphology  would  be  far  too  fictitious  to 
merit  use. 

For  each  channel  that  can  be  traced  (e.g.,  Figure  2),  I  will  transform  all  seismic  data  crossing  that 
channel  into  a  new  coordinate  system  consisting  of  length  along  and  distance  across  the  thalweg  of  the 
channel,  which  itself  must  be  interpolated  from  limited  sampling.  This  will  remove  the  sinuosity  of  the 
channel  from  the  problem,  and  allow  interpolation  by  conventional  geostatistical  means.  Transforming 
back  into  real  space  will  complete  the  full  interpolation  of  the  channel.  Merging  of  channels  can  most 
probably  be  accomplished  by  simple  averaging.  The  New  Jersey  seismic  data,  which  are  fully 
processed  and  in  hand,  will  serve  as  the  test  case  for  the  channel  interpolation  algorithm.  In  addition, 
the  “R”  and  “T”  horizons  will  be  interpolated,  using  methods  already  established,  to  provide  a 
complete  stratigraphic  model  of  the  area  of  importance  to  the  Geoclutter  program. 

WORK  COMPLETED 

Work  on  this  project  has  not  yet  commenced.  Work  will  begin  in  the  fall  and  be  completed  by  the 
summer  of  2001. 

RESULTS 

Results  have  not  yet  been  obtained  for  this  project. 

IMPACT/APPLICATIONS 

This  work  should  have  immediate  impact  in  the  Geoclutter  program  by  fully  enabling  acoustic 
reverberation  modeling  efforts  in  three  dimensions,  rather  than  constrained  in  2-dimensions  to  the 
orientation  of  available  seismic  profiles.  As  more  geophysical  data  are  collected  by  the  Geoclutter 
program,  these  products  can  be  incorporated  to  produce  a  refined  and  more  accurate  stratigraphic 
model  for  the  region. 

TRANSITIONS 

The  SimStrat  and  seafloor  conditional  simulation  programs  are  being  considered  by  the  Navy  for 
applications  in  minimum  data  density  analysis. 
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RELATED  PROJECTS 


The  geologic  and  geophysical  components  of  the  Geoclutter  program  are  a  direct  outgrowth  of  the 
ONR  STRATAFORM  program. 
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Figure  1.  Huntec  seismic  profile  on  the  New  Jersey  middle  shelf,  in  ~80  m  of 
water.  Horizon  interpretations  are  from  Duncan  et  al.,  2000. 
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"Channels"  Map  (2D  Huntec) 
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Figure  2:  Map  of  “ channels ”  horizon  on  the  New  Jersey  shelf  overlain  on  STRATAFORM  bathymetry.  Dashed  lines 
show  interpretation  of  the  drainage  system.  Note  complete  lack  of  correlation  between  bathymetry  and  channels 

morphology.  From  Duncan  et  al.  (2000). 


